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Abstract 

The relativistic generalization of the Brownian motion is discussed. We show that 
the transformation property of the noise term is determined by requiring for the 
equilibrium distribution function to be Lorentz invariant, such as the Jiittner dis- 
tribution function. It is shown that this requirement generates an entanglement 
O-ij between the force term and the noise so that the noise itself should not be a covari- 

r~| ' ant quantity. 
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1 Introduction 



Presently the relativistic hydrodynamics is found to be one of the important 
tool to describe the collective behaviors of relativistic heavy- ion collisions (jl|). 
However, there are several problems yet to be clarified when one wants to ex- 
tract quantitative information of the QGP formed in the collision (0). Some of 
them are associated with the difficulties of obtaining physical inputs such as 
equation of states, initial conditions or hadronization processes. Other ques- 
tions are associated with the basic principle of hydrodynamics itself, that is, 
the validity of the assumption of the local thermal equilibrium. Thus in order 
to extract conclusive information on the dynamics of relativistic heavy-ion 
collisions, we need to study these questions from various point of view. For 
example, to investigate more detailed behaviors than the simplest ideal fluid 
models and sudden freezeout mechanism, we need to develop mesoscopic the- 
ories like the Boltzmann equation, the Langevin equation and so on (0), to 
connect the hydrodynamic model to real observables. 

The relativistic Boltzmann equation approach has been applied by several au- 
thors to estimate the viscosity, heat conduction, etc (jl). Strictly speaking, the 
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predictions besed on the Boltzmann equation is reliable only in the Boltzmann- 
Grad limit and its applicability seems to be quite limited to extract the correct 
values of required observables in many realistic situations. 



As one possible alternative, we may consider the Langevin approach. This 
approach has been studied in low energy heavy-ion physics^) and also recently 
in relativistic cases particularly related to the behavior of heavy quarks (0). In 
the latter cases, the relativistic covariance of the Langevin approach deserves 
a careful attention. To the authors' knowledge, there are a very few studies 
on this aspect of the Langevin equation (0; 1(J 11; 12; 13; 14; 15; [l6l; [l7l). 
and the conclusions achieved in the former studies seem to raise some basic 
questions. 



In this paper, we would like to discuss the relativistic generalization of the 
Brownian motion. Actually the requirement of the relativistic covariance casts 
non-trivial questions because of the following reasons. The first point is the 
non-uniqueness of the discretization scheme of the Langevin equation due to 
the nature of multiplicative noise of a covariant Langevin equation for a rel- 
ativistic particle. Even if the noise is additive in a certain frame, it becomes 
multiplicative noise in other frame by the Lorentz boost. As is well known 
(and we will see later), when we have a multiplicative noise, the discretiza- 
tion scheme of the noise term is not unique and the solution depends on 
the descretization scheme. That is, the discretization scheme and the Lorentz 
transformation entangle each other. 

The second point is the transformation property of the noise. In general, the 
Langevin equation is composed by the force term and the noise term. In most 
of works, the force term and the noise term are assumed to be independent 
Lorentz vectors, and thus transform independently. In the continuum limit, 
this could be a natural assumption as far as the time evolution is defined in 
terms of the proper time of the particle. However, on integration, we have to 
define the Stieltjes integral of the noise term on a finite integration measure of 
the time. As a result, if the Lorentz covariance is required for the noise term, 
then the equilibrium distribution function depends on the integrating Lorentz 
frame and contradicts with the scalar nature of the equilibrium distribution 
function. 



In the following, we discuss the consistent relativistic generalization of Brow- 
nian motion. In Sec. 2, we first discuss the Brownian motion of a relativistic 
particle in the rest frame of a heat bath. We then analyze the corresponding 
Fokker-Planck equation by introducing the several integration schemes of the 
noise term. An explicit form of the distribution function at thermal equilib- 
rium is obtained. In Sec. 3, we show that how the Langevin equation obtained 
should transform under an arbitrary Lorentz boost of the system. This is 
done by requiring that the equilibrium distribution should be a scalar. The 
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result show that the Langevin equation of the relativistic Brownian motion 
has a non-trivial Lorentz transformation property. Summary and concluding 
remarks are given in Sec. 4. 



2 Relativistic Brownian motion in the rest frame of heat bath 



We consider the Brownian motion of a relativistic particle with mass m in the 
3+1 dimension. In the rest frame of a heat bath, we consider the following 
Langevin equation (in the units of c = 1), 



dt* p *' ^ ) 

% = -v(p°*)p* + y/2D(jP*)F(f), (2) 

where p°* = ^ (p*) 2 + m 2 . The parameter uijp*) and D{jP*) are assumed the 
Lorentz scalar functions and characterizes the relaxation of the momentum and 
the strength of the noise, respectively. Here the index * denotes the variables 
in the rest frame of the heat bath. The Gaussian white noise F(t) has the 
following correlation properties, 

<F(f*)>o = 0, (3) 
(F*(OF^*')) = <M(** -**')■ (4) 

The symbol {X) denotes the stochastic average of X in the rest frame of the 
heat bath (we refer to as RF). 

Now we replace the Langevin equation with the stochastic differential equation 
(SDE) 



dx* = ^dt\ (5) 
dp* = -u(p°*)p*dt* + ^/2D{p *)dw t * . (6) 
Here we used 

dwf = -w t *+dt* ~ w t . = F(t*)dt*, (7) 
with the help of the Wiener process w t * . The correlations are given by 
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(dw*,) = 0, 



(9) 



The last term of Eq. (jfJJ) contains a kind of Stieltjes integral of the stochastic 
variable. However, the definition of the required Stieltjes integral is not unique. 
Here we consider the three typical cases. 



^1) Ito interpretation (1181 ) 

In this case, the last term is interpreted as 



2£>(p°*)dw*. — ipD(pV*)dw* t *}j = pD( p v*(t*))(w* t * +dt 



(10) 



(2) Stratonovich-Fisk interpretation (I la ) 

In this case, the last term is interpreted as 



2D(p°*)dw t * — > 

J2D(p°*(t* + dt*)) + J2D(p°*(t*)) 



'11} 



By using the Ito formula (See appendix), we can show that the Stratonovich- 
Fisk SDE is equivalent to the following Ito SDE, 



dx.* = ^df 



:i2) 



dp* = -z/(p°*)p*dt* + yJD(f*)(d p .yjD(zP*))dt* + [^2D(p°*)dw* t ,}j. 

(13) 



(3) Hanggi-Klimontovich (UJ; |20|; M) 

In this case, the last term is interpreted as 



2D(p°*)dw t * 



2D(p°*)dw* t *] HK = yj2D(pP*(t* + dt*))(w* t , +dt , - w**) 

(14) 



By using the Ito formula, we can show that the Hanggi-Klimontovich 
SDE is equivalent to the following Ito SDE, 



P 

dx = —z-dt 



P 



o* 



(15) 



dp* = -is(p v *)p*dt* + 2^/D(p *)(d p ^D(p°*))dt* + [^2D(p°*)dw*,]i. 

(16) 
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We are interested in the equilibrium distribution function described by using 
these SDEs. For this purpose, we introduce the probability density 

p(x, P ,t) = (^)(x-x*(r))5( 3 )( P - P *(r))). (17) 



Then, the time evolution of p(x, p, t) is given by the Fokker-Planck equation 
is 



(W(x*, p*, O = - E ^p(x*, pV*)) + E Kp°V p*, O) 

i P i 

+ E J D 1 - a ( P *)^D a (p°*)p(x*, P *, e). (is) 



The parameter a denotes the different discritization scheme, a = 0, 1/2 and 
1 represents the Hanggi-Klimontovich, Stratonovich-Fisk and Ito scheme, re- 
spectively. 

The equilibrium distribution of the Fokker-Planck equation is 



p st (p*) oc exp 



D(q 



(19) 



) 



where the path C of the integral in the momentum space is arbitrary. This 
equilibrium distribution should be a Lorentz scalar by definition Eq. (fl7l) . As 
we will see later, we can define the transformation property of the noise by 
using this fact. 



3 Relativistic Brownian motion in general frame 



We consider the reference frame which is moving with the velocity V with 
respect to the rest frame of the heat bath (we refer to as simply MF-moving 
frame). The four-momentum djf in this frame is then given by the Lorentz 
transformation of dp* 11 as 



dp^ = A(V)dp*^ 



( 



7 (V) /3(V> r 7 (V) 
/3(V> 7 (V) 7 (V)P||+Q 




(20) 



5 



Here, Pu = nn T and Q± = i — P\\ with n = V/|V|. By using the on mass-shell 
condition 



dp° = ~^dp 



p° 



we get the SDE in the MF as 

dpl = _ H^PMV) |^ 0(pJ _ ^(y^ipO) + f3(V)(p 2 n l - pvtf)} dt + [Bdw* t *)\ (21) 



where pv = n T p and 



B = y/2D(u?pJ { P °P|| + 7 (V) (p° - (3(V)p v + 0(V)(n ■ p T )) Q ± } . 

(22) 



The last term in Eq. (12 lj) should be interpreted appropriately according to the 
integration scheme, namely the Ito, Stratonovich-Fisk or Hanggi-Klimontovich 
cases. 

The last term is not yet transformed because it contains the noise in RF dw t * . 
First of all, we assume that, even after the Lorentz boost, the stochastic part 
of the Brownian motion still preserves the property of the Gaussain white 
noise, which is defined by 



(dw t ) v = 0, (23) 
(dw i t dwij v = dt5 ij 5i m , (24) 

Here the simbol (X) v denotes the stochastic average of X in the MF. 

Then, at first sight, as was assumed in previous works (J7J; 

[2 QBE El), 

we may suppose that the noise term is isolatedly Lorentz 
covariant and we could write 



(dw t )v=(dw*) = 0, (25) 
7- 1 (p)(rfw^w|> = 7(p*)(rfw|dwf } = drSuSu, (26) 



where 



r\p) = (A(V)A(-p*)) 00 = 7 (V) 7 (P*)(1 - P(V)^). (27) 
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If this is the case, we could define the transformation property of the noise as 
follows; 



dw* t , = ^dw t = ^/ 7 (p)7-i(p*)dw, = i l/2 iy)^ V pi )Pv dw t . (28) 

However, as we will see in below, this transformation rule (|28|) does not give 
the correct properties of the Langevin equation. The reason for this is that the 
Stieltjes integral associated with the noise term is defined on the time interval 
dt, so that dw*» is non-local in the time t. Thus the Lorentz transformation 
entangles with the integration scheme in the order of dt. Then the noise term 
itself is not covariant but constitutes a Lorentz vector together with the force 
term. That is, the force part and the stochastic part can be mixed. Then, the 
first order correlation calculated in the MF of the noise term in the rest frame 
of the heat bath does not vanish, 



(dw* r ) v ^0. (29) 



Thus we should consider, instead of Eq. fl28|) the following transformation prop- 
erty of the noise, 



dt* 

dw% = \ ——dw t + C u dt 
t » dt p 



= inV) f-^ V)PV d* t + C,d t . (30) 

Here C p dt term acts as the force term in the MF, which should be separated 
from the pure stochastic part dw t . By using this definition and using the Ito 
formula, the Langevin equation is given by 



dp* = {pV - /3(V)ny) + (3(V)(p 2 n l - Pv p l )}dt 

+ (1 - a) ]T W k dl,B ik dt + {BC p fdt + [Bdw t ]i, (31) 

jk 



where 



B = pg3B B (32) 

and the last term should be calculated according to the Ito scheme. It should 
be noted that, to obtain a consistent result, the Ito formula to convert the 
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SDE from one scheme to the other and the replacement of the noise by using 
Eq. (|30|) must be used only after the transformation of the noise, Eq. ( 130]) . 



Let us calculate the equilibrium distribution function in the MF. It is the 
solution of the equation, 



-A* - (1 - a) B jk dlB ik + -<9 P (BB T ) ij ' 

jk 



Pst(p) = 0, (33) 



where 



A = -^pl | p o (pJ _ p(y) n i p 0) + f3(v)(p 2 n l - p vV 1 )} + [BC p }\ (34) 

The logarithmic derivative of the equilibrium distribution with respect to the 
momentum is given by 

(dlPst)/Pst = 2(BB r )r 1 (A* + (1 - a)(d l p W k )B lk - \d k p (BB T f\ .(35) 



If p st is a scalar function, it should coincide with the logarithmic derivative 
of the same equilibrium distribution function obtained in the RF. From Eq. 
( TT91) . we have 

/ pn T Y k f is(p *) D'(p°*)v* 



where D'(p°) = dD(p°)/dp°. To compare Eq. (l35l) with Eq. (l36|) . it is convenient 
to decompose the vector C p into two directions, 



C p = C£(p-p v n) + Cjln. (37) 
Then 



C; = -(l-a) 



D(n^)/j(\/)7(K)^-/3(\/)p 
2 (p ) 2 p°-/3(V>v : 



l D{u»p„) 1 /?(l / )m 2 



2 (/) 2 p°-/3(% 

(^2 



+(l-2a) \/ 7- n ,„ (p 



+(1 



l D(vMp„)P(V) 7 (V) jfp 2 \ fo ^ l p v -(3(V)p° \ 

(39) 



where d denotes the dimension of the space part. One can see that the form of 
C p depends only on the parameter of the strength of the noise D(p *), although 
the original Langevin equation contains two parameters, u(p *) and D(p *). 
We further notice that the C p does not vanish by changing D(p *). That is, in 
order to keep the scalar property of the equilibrium distribution function, C p 
cannot be null. It means that we cannot use the commonly assumed Lorentz 
covariant noise Eq. (128]) . 



When the equilibrium distribution function is given by the Juttner function, 

p st = Canst, x e-^", (40) 
then the parameters of the SDE satisfy the following relation, 

f(u"p**) = -^(/^KPm) - aD'(u^)), (41) 
u^p^ 



which is the generalized Einstein's dissipation-fluctuation relation of the rela- 
tivistic Brownian motion. 



4 Summary and Concluding remarks 



In this work, we discussed the generalization of the Brownian motion of a 
relativistic particle. First, we prepared a SDE in the rest frame of a heat 
bath. We then derived the equilibrium distribution function by calculating 
the corresponding Fokker-Planck equation. Next, by carrying out the Lorentz 
boost explicitly requiring the on-mass-shell condition, we obtained the boosted 
SDE in the moving frame with respect to the heat bath. The covariance of the 
SDE requires that the noise is essentially multiplicative in general frame. We 
also found that the Lorentz boost induces a non-trivial entanglement between 
the force term and the noise term. We demonstrated that the commonly used 
Lorentz invariant noise does not lead to an invariant equilibrium distribution 
in the present formulation of the relativistic Brownian motion. 



The equation discussed here is different from those proposed in (Il4l;ll6l;ll7l). In 
the RF, our equation is equivalent to their equations. However, they assume 
the existence of the Lorentz covariant noise term, and their equations in the 
MF are different from ours. To our best knowledge, the consistency of the 
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equilibrium distribution functions in the RF and the MF is not discussed in 
their approaches. 



The present study will be useful to investigate the Brownian motion of a 
relativistic particle in a relativistically expanding heat bath such as jet propa- 
gation in an expanding quark-gluon plasma formed in ultra-relativistic heavy 
ion collisions. The further studies are in progress. 
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A Ito formula 



Let us consider an arbitrary function /(x) and the evolution of x is given by 
the SDE, 



rfx = Adt + [Bdw]/. (A.l) 
Then, the variation of /(x) is 

d/(x) = EA'fik/Cx) + ^ ^PB^M/WK + E W/(x)^]/.(A.2) 

i ij ij 



This is called the Ito formula (I la ). By using the Ito formula, we obtain, 



[Bdw] SF = [Bdw}\ + \Y. H^dj^dt, (A.3) 

jk 

[Bdw} j HK = [Bdw]) + £ B jk djB ik dt. (A.4) 

jk 

Thus we can conclude as follows. When we have the Stratonovich-Fisk SDE, 



rfx = Adt + [Bdw] SF , (A. 5) 

this is equivalent to the Ito SDE, 
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dx* = {A i + - B jk djB ik }dt + [Bdw]/. (A.6) 

jk 

When we have the Hanggi-Klimontovich SDE, 



dx = Adt + [Bdw] HK , (A. 7) 

this is equivalent to the Ito SDE 



dx* = {A' + J2 B j %B ik }dt + [Bdw]/. (A.8) 

jk 
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